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The meson masses and mixing angles have been calculated for the scalar and pseudoscalar sector 
in the framework of the generalized 2+1 flavor Polyakov loop augmented quark meson linear sigma 
model. We have given the results for two different forms of the effective Polyakov loop potential. 
The comparison of results with the existing calculations in the bare 2 + 1 quark meson linear sigma 
| model, shows that the restoration of chiral symmetry becomes sharper due to the influence of the 

Polyakov loop potential. We find that inclusion of the Polyakov loop in quark meson linear sigma 
■ model together with the presence of axial anomaly, triggers an early and significant melting of the 

strange condensate. We have examined how the inclusion of the Polyakov loop qualitatively and 
quantitatively affects the convergence in the masses of the chiral partners in pseudoscalar (n, rj, rj' , 
qq ■ K) and scalar (<r, eio, /o,«0 meson nonets as the temperature is varied on the reduced temperature 

scale. The role of Ua(1) anomaly in determining the isoscalar masses and mixing angles for the 
pseudoscalar (77 and 77') and scalar (a and /o) meson complex, has also been investigated in the 
Polyakov quark meson linear sigma model. The interplay of chiral symmetry restoration effects and 
i-C ' the setting up of Ua(1) restoration trend has been discussed and analyzed in the framework of the 

presented model calculations. 
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> ' I. INTRODUCTION 

^ . 

, The present theoretical understanding of strong interaction physics indicates that normal hadronic matter undergoes 
^| ' a phase transition, where the individual hadrons dissolve into their constituents and produce a collective form of matter 
known as the Quark Gluon Plasma (QGP) under the extreme conditions of high temperature and/or density PHI]- 
Rclativistic heavy ion collision experiments at RHIC (BNL), LHC (CERN) and the future CBM experiments at the 
' FAIR facility (GSI-Darmstadt) aim to create and study such a collective state of matter. 

Study of the different aspects of this phase transition, is a tough and challenging task because Quantum 
' Chromodynamics (QCD), the theory of strong interaction, becomes nonperturbative in the low energy limit. In 
the zero quark mass limit, chiral condensate works as an order parameter for the spontaneous breakdown of the chiral 
symmetry in the low energy hadronic vacuum of the QCD. For the infinitely heavy quarks, in the pure gauge SU C (3) 
QCD, the Z(3) (center symmetry of the QCD color gauge group) symmetry, which is the symmetry of hadronic 
vacuum, gets spontaneously broken in the high temperature/density regime of QGP. Here the expectation value of 
the Wilson line (Polyakov loop) is related to the free energy of a static color charge, hence it serves as the order 
parameter of the confinement-deconfinement phase transition Q ■ Even though the center symmetry is always broken 
with the inclusion of dynamical quarks in the system, one can regard the Polyakov loop as an approximate order 
parameter because it is a good indicator of the confinement-deconfinement transition Q . 

We get important information and insights from the lattice QCD calculations (see e.g. 043) regarding various 
aspects of the transition, like the restoration of chiral symmetry in QCD, order of the confinement-deconfinement 
phase transition, richness of the QCD phase structure and mapping of the phase diagram. Since lattice calculations 
are technically involved and various issues are not conclusively settled within the lattice community, one resorts to the 
calculations within the ambit of phenomcnological models developed in terms of effective degrees of freedom. These 
models serve to complement the lattice simulations and give much needed insight about the regions of phase diagram 
inaccessible to lattice simulations. A lot of current effective model building activity, is centered around combining the 
features of spontaneous breakdown of both chiral symmetry as well as the center Z(3) symmetry of QCD in one single 
model (see for example [l5l - l29j ). In these models chiral condensate and Polyakov loop are simultaneously coupled to 
the quark degrees of freedom. 
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In order to calculate the properties of mesons in hot and dense medium, several investigations have been done 
in the two and three flavor Nambu-Jona-Lasinio (NJL), Polyakov-Nambu-Jona-Lasinio (PNJL) models (e.g. j30l - l33j ) 
and also in the SU(2) version of linear sigma model (e.g. |34l - l36l |). Since chiral symmetry restoration is signaled by 
parity doubling, these studies look for the patterns of emerging convergence in the masses of the chiral partners in 
pseudoscalar (jr, 77, 77', K) and scalar mesons (a, ao, fo, k). It is a common knowledge that the basic QCD Lagrangian 
has the global SUl+r(3) x SUl-r(3) x Ua{^) symmetry. Different patterns of spontaneous as well as explicit breaking 
of SUv(3) x SUa(3), have been discussed by Lenaghan et al. [37} in the ambit of SU(3) linear sigma model. Schaefer 
et al. enlarged the linear sigma model with the inclusion of quarks [381 ] and then they studied in the 2+1 flavor 
breaking scenario, the consequences of SU(3) chiral symmetry restoration for scalar and pseudoscalar meson masses 
and mixing angles, in the presence as well as the absence of 1/4(1) axial symmetry, as the temperature is increased 
through the phase transition temperature. The Ua{^) axial symmetry does not exist at the quantum level and as 
shown by 't Hooft [39) . it gets explicitly broken to ZA(Nf) by the instanton effects. The J7a(1) anomaly does not let 
the 77' meson remain massless Goldstone boson in the chiral limit by giving it a mass of about 1 GeV. This happens 
due to flavor mixing, a phenomenon that lifts the degeneracy between the 7r and 7/ which otherwise would have been 
degenerate with 7r in £/(3) even if the explicit chiral symmetry breaking is present. There is large violation in the 
Okubo-Zweig-Iizuka rule for both pseudoscalar and scalar mesons and ideal mixing is not achieved because of strong 
flavor mixing between nonstrange and strange flavor components of the mesons [31| . Hence J7a(1) restoration will 
have important observable effects on scalar and pseudoscalar meson masses as well as the mixing angles. 

In a three flavor PNJL model calculation, Costa and collaborators [3l[ have discussed in detail how the inclusion 
of Polyakov loop in the NJL model, affects the results of meson mass and mixing angle calculations. However in an 
earlier paper, they have pointed out that the description of the 77' in the NJL model has some problem [32|]. The 
NJL model does not confine and the meson degrees of freedom are generated in the model by some prescription. 
The polarization function for the meson gets an imaginary part above the qq threshold, hence 77' becomes unbound 
completely in the model soon after the temperature is raised from zero. Thus 77' in the NJL model is not a well defined 
quantity [4(| ■ Schaefer et al. [H, |4l| have also made an elaborate study of meson masses and mixing angles with and 
without Ua{^) axial anomaly in the 2+1 flavor quark meson linear sigma model where the mesons are included in 
the Lagrangian from the very outset and the J7a(1) breaking 't Hooft coupling term is constant. The behavior of the 
scalar and pseudoscalar mixing angles in their calculation is opposite to what has been reported in the calculation 



by Costa et al. [311 ] . It is worthwhile and important to investigate the influence of Polyakov loop on meson mass 



and mixing angle calculations in scalar and pseudoscalar sector, in the framework of generalized 2+1 flavor quark 
meson linear sigma model enlarged with the inclusion of the Polyakov loop [i^l - liij . Since we are lacking in the 
experimental information on the behavior of mass and mixing angle observables in the medium, a comparative study 
of these quantities in different models and circumstances becomes all the more desirable. We will be investigating 
how the inclusion of Polyakov loop, qualitatively and quantitatively affects the convergence of the masses of chiral 
partners, when the parity doubling takes place as the temperature is increased through T c and the partial restoration 
of chiral symmetry is achieved. We will also be studying the effect of Polyakov loop on the interplay of SUa(3) chiral 
symmetry and Ua{X) symmetry restoration. 

The arrangement of this paper is as follows. In Sec|H]we have given the formulation of the model. The description 
of grand potential in the mean field approach has been presented in Sec. IIIII We have derived the modification of 
meson masses due to the qq contribution in the presence of Polyakov loop in Sec lIVI where the formulae for meson 
masses and mixing angles have been discussed. In Sec|Vl we will be discussing the numerical results and plots for 
understanding and analyzing the effect of Polyakov loop on chiral symmetry restoration. Discussion and summary is 
presented in the last Sec lVIl 



II. MODEL FORMULATION 



We will be working in the generalized three flavor quark meson linear sigma model which has been combined with 
the Polyakov loop potential [4^h44| . In this model, quarks coming in three flavor are coupled to the SUv{3) x SUa(3) 
symmetric mesonic fields together with spatially constant temporal gauge field represented by Polyakov loop potential. 
Polyakov loop field $(x) is defined as the thermal expectation value of color trace of Wilson loop in temporal direction 



$ = — Tr c i, 



+ 



1 



Tr r tf 



(1) 



where L(x) is a matrix in the fundamental representation of the SU C (3) color gauge group. 



L(x) = Vexp 



cLtAq (x, t) 



(2) 
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Here V is path ordering, Ao is the temporal component of Euclidean vector field and j3 = T 1 |4| . 

The model Lagrangian is written in terms of quarks, mesons, couplings and Polyakov loop potential U ($, $*, T). 

C pqms = Cqms-U($,®*,T) (3) 
where the Lagrangian in quark meson linear sigma model 

C-QMS = <Tf (i7 M - D M - 9 T a + H5^a) ) Qf + C m (4) 

The coupling of quarks with the uniform temporal background gauge field is effected by the following replacement 
Dfj, = dp iAp and = S^Aq (Polyakov gauge), where A^ = g s A a ^X a /2. g s is the SU C (S) gauge coupling. A a are 
Gell-Mann matrices in the color space, a runs from 1 • • • 8. qf = (u,d, s) T denotes the quarks coming in three flavors 
and three colors, g is the flavor blind Yukawa coupling that couples the three flavor of quarks with nine mesons in 
the scalar (a a , J p = + ) and pseudoscalar (ir a , J p = 0~) sectors. 

The quarks have no intrinsic mass but become massive after spontaneous chiral symmetry breaking because of 
nonvanishing vacuum expectation value of the chiral condensate. The mesonic part of the Lagrangian has the following 
form 

C m = Tr (d^d^M) ~ m 2 Tr(M t M) - Ai [Tr^M)] 2 
— A 2 Tr (M f M) 2 + c[det(M) + det(M^)} 

+Tr [H(M + M + )] . (5) 
The chiral field M is a 3 x 3 complex matrix comprising of the nine scalars a a and the nine pseudoscalar 7r a mesons. 

M = T a £ a = T a (a a + in a ) (6) 
Here T a represent 9 generators of U(3) with T a = 4f . a = 0, 1 . . . 8. A a are standard Gell-Mann matrices with 



A = yj 3 1. The generators follow U(3) algebra [T ai T b ] = if a bcT c and {T a ,T b } = d abc T c where f abc and d abc are 

standard antisymmetric and symmetric structure constants respectively with f ab o = and d a bo = yf 1 <^ab and 

matrices are normalized as Tr (T a T b ) — 

The SUl(3) x SUr(3) chiral symmetry is explicitly broken by the explicit symmetry breaking term 

H = T a h a (7) 

Here H is a 3 x 3 matrix with nine external parameters. The £ field picks up the nonzero vacuum expectation value, £ 
due to the spontaneous breakdown of the chiral symmetry. Since £ must have the quantum numbers of the vacuum, 
explicit breakdown of the chiral symmetry is only possible with three nonzero parameters ho, /13 and h%- We are 
neglecting isospin symmetry breaking hence we choose ho, /is 7^ 0. This leads to the 2 + 1 flavor symmetry breaking 
scenario with nonzero condensates ao and as- 

Apart from ho and hg, the other parameters in the model are five in number. These are the squared tree- level mass 
of the meson fields m 2 , quartic coupling constants Ai and A2, a Yukawa coupling g and a cubic coupling constant c 
which models the Ua{^) axial anomaly of the QCD vacuum. 

Since it is broken by the quantum effects, the J7a(1) axial which otherwise is a symmetry of the classical Lagrangian, 
becomes anomalous J45| and gives large mass to 77' meson (m n i = 940 MeV). In the absence of Ua(1) anomaly, r/' 
meson would have been the ninth pseudoscalar Goldstone boson, resulting due to the spontaneous break down of the 
chiral Ua{3) symmetry. The entire pseudoscalar nonet corresponding to spontaneously broken 6^(3), would consist 
of the three tt, four K, rj and rj mesons, which are the massless pure Goldstone modes when H = and they become 
pseudo Goldstone modes after acquiring finite mass due to nonzero H in different symmetry breaking scenarios. The 
particles coming from octet (ao, fo, k) and singlet (a) representations of SUy(3) group, constitute scalar nonet (cr, 
ao, fo, k). In order to study the chiral symmetry restoration at high temperatures, we will be investigating the trend 
of convergence in the masses of chiral partners occurring in pseudoscalar (n, r\, 77', K) and scalar (cr, ao, fo, «) nonets, 
in the 2 + 1 flavor symmetry breaking scenario. 



A. Choice of Potentials for the Polyakov Loop 



The effective potential U ($,$*, T) is constructed such that it reproduces thermodynamics of pure glue theory 
on the lattice for temperatures upto about twice the deconfinement phase transition temperature. At much higher 
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TABLE I: The squared masses of scalar and pseudoscalar mesons appear in nonstrange-strange basis. In the following table x 
denotes u x and y denotes <r y . The masses of nonstrange <tns, strange as, nonstrange rjNS and strange rjs mesons are given in 
the last two rows. 



temperatures, the transverse gluons become effective degrees of freedom, hence the construction of effective potential 
in terms of the Polyakov loop potential is not reliable [3, H(| ■ 

At low temperatures, the effective potential U ($,$*, T) has only one minimum at $ = in the confined phase. 
Above the critical temperature for deconfinement transition, $ = becomes metastable local minimum and now, the 
effective potential has three degenerate global minima at $ ^ due to the spontaneous breakdown of the Z(3) center 
symmetry. 

In this work, we use the following two choices of the effective potential. The first choice is based on the polynomial 
expansion in terms of Polyakov loop order parameter $ and is given [l8| as 



T 4 4 V J 6 

+ ^(m 2 + \^\ 2 ) 2 (8) 

The second term that is the sum of $ 3 and <I>* 3 terms, causes the three degenerate vacua above the deconfinement 
phase transition. The potential parameters are adjusted according to the pure gauge lattice data such that the 
equation of state and Polyakov loop expectation values are reproduced. The temperature dependent coefficient 62 (T) 
governs the confincmcnt-dcconfinement phase transition and is given by 

b 2 (T) =a + ai +a 2 + ° 3 ( Y 

The other parameters have the following value 

a = 6.75, ai = -1.95, a 2 = 2.625 , 
a 3 = -7.44, 63 = 0.75 , 64 = 7.5 . 

The other choice of effective potential as given in Ref. , has the logarithmic form. The results produced by this 
potential are known to be fitted well to lattice results. 



W log ($,$*,T) a(T) 



$*$ + 6(T) ln[l-6$*$ 



T 4 2 

+4($* 3 + $ 3 ) -3($*$) 2 ] (9) 



where the temperature dependent coefficients arc as follow 

a(T) = ao + ai (|)+ a2 (|) 2 b(T) = b 3 (^ 3 

The critical temperature for deconfinement phase transition To = 270 McV is fixed for pure gauge sector. The 
parameters of Eq. © are 

00=3.51, ai = -2.47, 
a 2 = 15.2 , b 3 = -1.75 
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Both effective potential fits reproduce equally well the equation of state and the Polyakov loop expectation value. 





C[MeV] 


m 2 [MeV 2 ] 
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(342. 52) 2 
-(189.85) 2 
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TABLE II: parameters for m a — 600 MeV with and without £/a(1) axial anomaly term. 



III. GRAND POTENTIAL IN THE MEAN-FIELD APPROACH 



The thermodynamics of changing numbers of particles and antiparticles is governed by grand canonical partition 
function. We are considering a spatially uniform system in thermal equilibrium at finite temperature T and quark 
chemical potential /Lt/(/ = u,d,s). The partition function is written as the path integral over quark/ antiquark and 
meson fields I3a| 



Z = Trcxph/?(H- 



f—u,d,s 



= J \{V(T a Vir a J VqVq exp - jf dr J d' 



f=u,d,s 



(10) 



where V is the three dimensional volume of the system, and ft = y. For three quark flavors, in general, the three 
quark chemical potential are different. In this work, we assume that S'L r y(2) symmetry is preserved and neglect the 
small difference in masses of u and d quarks. Thus the quark chemical potential for u and d quarks become equal 
fJ'x = li-u = \id- The strange quark chemical potential is fj, y = fi s . Further we consider symmetric quark matter and 
net baryon number to be zero. 

The partition function for the SU(3) version of the linear sigma model with or without quarks can be evaluated 
by the more advanced many-body resummation techniques such as the self consistent Hartree approximation in the 
Cornwall, Jackiw and Tomboulis [37l |46| formalism or the so called optimized perturbation theory |47j in its improved 
version j48l450| . However the predictive power of these methods depends on how they are implemented in different 
approximation schemes. 

In the simple mean field approximation, one does not encounter various problems of the more advanced many-body 
resummation techniques. In our work, the partition function has been evaluated in the mean-field approximation 
We replace meson field by their expectation values ($) = Tqctq + T%a% and neglect both thermal as well as 
quantum fluctuations of meson fields while quar ks and antiquarks are retained as quantum field. Now following the 
standard procedure as given in Refs. (TBI . Il8l. l5ll. l52j one can obtain the expression of grand potential as sum of pure 
gauge field contribution IA ($,$*, T), meson contribution and quark/antiquark contribution evaluated in the presence 
of the Polyakov loop, 



Q(T,/i) = -^^ = U(a ,a s )+U (<&,$*, T) 



V 



(11) 



In order to study 2 + 1 flavor case, one performs following basis transformation of condensates and external fields 
from original singlet octet (0, 8) basis to nonstrange strange basis (x, y). 



(12) 
(13) 
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Similar expressions exist for writing the external fields (h x , h y ) in terms of (ho, hg). Thus the nonstrange and 
strange quark/antiquark decouple and the quark masses become 



9 Y , 



V5 



(14) 



Quarks become massive in symmetry broken phase because of non zero vacuum expectation values of the condensates. 
The mesonic potential in the nonstrangc-strangc basis reads, 



U(a x ,a y ) 



h x cr x 



kyCTy 



+ -(2A 1+ 2A 2 )^ 



2^2 



(15) 



The chiral part of the Polyakov loop augented quark meson linear sigma (PQMS) model has the six input parameters 
and therefore require six known quantities as input. In general m^, tuk, the pion and kaon decay constant f w , fx, 
mass square of 77, rf and m a are used to fix these parameters. The parameters are fitted such that in vacuum the 
model produces observed pion mass 138 MeV. In the present work we are using the set of parameters for sigma mass 
m a = 600 MeV. The parameters used in this work, taken from |38[ , are shown in Table III Al 

Finally the quark/antiquark Polyakov loop contribution reads, 



-2T £ 

f—u,d,s 



d 3 p 



(27r) a 



ig f + ln -9/ 



(16) 



We define and after taking trace over color space 
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1 + 3<S>e~ E t ' T + 3$*e- 2 ^ /T + e^ E V T 



(17) 



9 f 



1 + Z<5>*e~ E } /T + Z<5>er 2E f /T + e" 3E J /T 



(18) 



Here we use the notation = Ef =F A* an d Ef i s the flavor dependent single particle energy of quark/antiquark. 



(19) 



m/ is flavor dependent quark mass and is function of condensates 00 and as- 
One can very easily notice from equations (jTTJ) and (TT5|) that the role of quarks and antiquarks as well as that of 
the Polyakov loop and its conjugate can be interchanged by the transformation fi — > —fx. Confinement is the very 
interesting feature of the QCD and the PQMS model describes this behavior qualitatively. The Polyakov loop is order 
parameter for confinement-deconfinement phase transition. In the confined phase $ = 0. It can also be noticed from 
the grand potential that one and two quark state contributions are vanishing. Only the three quark states contribute. 
In this way the PQMS model qualitatively mimics confinement of quark/antiquark within three quark color singlet 
states [22! ]. 

One can get the quark condensates a Xi a y and the Polyakov loop expectation values $, $* by searching the global 
minima of the grand potential for a given value of temperature T and chemical potential [i. 



da x 



on 

dcr v 



on 



on 



<9$* 



(20) 



In this work we are always considering the fi = case. 
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IV. MESON MASSES AND MIXING ANGLES 



The curvature of grand potential Eq. (fTTj) at the global minimum determines scalar and pseudoscalar meson masses. 



l a,ab 



Ct, I 



(21) 



where subscript a = s, p; s stands for scalar and p stands for pseudoscalar meson and a, b = ■ • • 8. We note 
that the Polyakov loop decouples from the mesonic sector at T=0 and the meson masses do not receive contribution 
from quark/ ant iquark in vacuum and hence meson masses are governed by mesonic potential only. The mesonic 
contribution to the meson masses is summarized in Table [TTJ The diagonalization of (0 - 8) component of mass matrix 
gives masses of a and /o mesons in scalar sector and masses of rj and rj in pseudoscalar sector. The scalar mixing 
angle 9 S and pseudoscalar mixing angle 9 P are given by, 



tan 28 r , 



2m l,os 



"a,00 



'a, 88 



(22) 



Here a stands for scalar and pseudoscalar field. The detail expressions for masses and mixing angles are given in 
Ref. (37l [lH . The meson masses are further modified in medium at finite temperature by the quark contributions in 
the grand potential. In order to calculate the second derivative Eq. (|2Tj) for evaluating the quark contribution in the 
presence of the Polyakov loop potential, the complete dependence of all scalar and pseudoscalar meson fields Eq.© 
has to be taken into account. We have to diagonalize the resulting quark mass matrix. The expression for the meson 
mass modification due to quark contribution at finite temperature in QMS model, has been evaluated by Schaefer et 
al. and is given as 



L a.ab 



(a++aj) 



L f,ab 



■■ E 

f=m,y 



d 3 p 1 



IE) 



f 



bj) 



m /,a m /,fe 

2E f T 



(23) 



dm 2 j I d£, a ,a is the first derivative and m 



with respect to meson fields £ Q! b. 
quark/antiquark occupation numbers; given as 



2 a /d£a,b is the second derivative of squared quark mass 
The number of internal quark degrees of freedom, v c = 2N C = 6. Here are 



f.ab = <>">). 



(24) 





2 2 / 
m x,a m x,b/ 


9 A m 2 x,ab/g 2 


m 2 y , a rnl, b /g A 


2 / 2 

m v ,ab/g 


To 


(To 


12 

3 x 


2 
3 


1 rfi 
3 a y 


1 

3 


<J\ 


(Tl 


1 „2 


1 








"\ 


<7-l 










y/2o- x +2a y 




y ^l-o'i 


<JX 
(TO 


(T* 
(T* 


1 „ 2 

6 x 


1 

3 

vl 
3 


_^2 CT 2 

3 V 


2 
3 

V2 
3 


7T0 


7To 





2 
3 





1 

3 


7T1 


TTl 





1 








714 


7T4 





cr x — y/2cr y 





■\/2o"a; — 2(Ty 






7Ts 


7T8 





1 

3 





2 
3 


7T0 


7T8 





3 





v/2 
3 



TABLE III: First and second derivative of squared quark mass in nonstrange-strange basis with respect to meson fields are 
evaluated at minimum. Sum over two light flavors, denoted by symbol x, are in third and fourth columns. The last two columns 
have only strange quark mass flavor denoted by the symbol y. 
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and the notations bj~a^ — (flj ) 2 stand for particle (+) and antiparticle ( — ) in quark meson linear sigma model 
without inclusion of the Polyakov loop. 

The expression of mass modification due to quark contribution at finite temperature, will change in the presence of 
the Polyakov loop. We are obtaining the following formula for the mass modification that results on account of quark 
contribution in the PQMS model 



^ m l,ab 



d 2 n qq (T,fi) 



(Aj + Aj)[mj iab 



3E 

f=x,y 



d 3 p 1 
(2^)3 E~ s 



m la m lb 



IE) 



(B++BJ) 



m ),a m ),b 

2E f T 



(25) 



The notations and By have the following definitions 
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4 



9f 



(26) 



$*e~ E f /T + 2<f>e~ 2E f /T + e~ 3E f /T 



9 f 



and B^ = 3(A^) 2 — C^, where we again define 



(27) 



<7 + = 



(28) 



c 7 = 



$* e ~ E f l T + 4$e _2 - E / /T + 3e~ 3E f /T 
9] 



(29) 



The squared quark mass derivatives evaluated at minimum which were originally derived in Ref . [38j , are collected 
in Table HTT1 The inclusion of Polyakov loop in QMS model does not make any change in these equations. 



V. EFFECT OF THE POLYAKOV LOOP ON THE RESTORATION OF CHIRAL SYMMETRY 



We are presenting the result of our calculation for estimating the effect of the Polyakov loop potential on the 
restoration of chiral symmetry when it is included in the 2 + 1 flavor quark meson linear sigma model at finite 
temperature and zero chemical potential with and without axial Ua(X) breaking. We have considered the two different 
ansatzs for the Polyakov loop potential namely the polynomial potential and log arithmic potential and compared the 
results with the existing calculations in the quark meson linear sigma model [38|. The interplay of the effect of Ua(X) 
axial restoration and chiral symmetry restoration in the presence of the Polyakov loop potential has been shown 
through the temperature variation of strange, nonstrange chiral condensates, meson masses and mixing angles. The 
Ua{]~) axial breaking term has been kept constant throughout the investigation. The value of Yukawa coupling g has 
been fixed from the nonstrange constituent quark mass m q = 300 MeV and is equal to 6.5. This predicts the strange 
quark mass m s ^ 433 MeV. 





QMS 


PQMS:pol 


PQMS:log 


T c x (MeV) 


146 


204 


206 


Tf (MeV) 


248 


262 


274 


Tf (MeV) 




204 


206 



TABLE IV: The characteristic temperature (pseudocritical temperature) for the chiral transition in the nonstrange sector T*, 
strange sector T s x and confinement-deconfinement transition T*, in the QMS, PQMS:log and PQMS:pol models. 
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FIG. 1: The variation of nonstrange a x , strange a y condensates with respect to the relative temperature scale (T/T*) at zero 
chemical potential (/i = 0) in the QMS model and PQMS models with polynomial and logarithmic potentials for the Polyakov 
loop is shown. The lines with continuous dots represent the variation in the QMS model, while the dashed-dotted lines show 
the variation in PQMS:pol model and the solid lines are the variations in the PQMS dog model. The line with the big solid dots 
shows the a y variation in the PQMSdog model, while the dark dashed line shows the pure QMS model results when anomaly 
is absent, i.e. c = 0. The expectation value of the Polyakov loop ($), in PQMS:pol and PQMSdog model is shown in the right 
plots. 



The solutions of the coupled gap equations, Eq. (|20p determine the nature of chiral and dcconfincmcnt phase 
transition through the temperature and chemical potential dependence of nonstrange and strange condensates (a x 
and a y ) and the expectation value of the Polyakov loop ((<&) and (*&*))• The temperature variation of a x , a v and 
($) in mean-held approximation, at zero chemical potential in the PQMS models with the polynomial Polyakov 
loop potential Eq.® (PQMS:pol) and the logarithmic Polyakov loop potential Eq.© (PQMSdog) is shown in FigQ] 
We have also plotted the strange and nonstrange condensate in QMS model to compare and investigate the effect 
of the Polyakov loop potential inclusion on chiral symmetry restoration trend reflected through masses and mixing 
angles of mcsonic excitations. The characteristic temperatures (pseudocritical temperature) for the confinement - 
deconfincment transition T*, the chiral transition in the nonstrange sector T* and strange sector are defined 
through the inflection point of ($), a x , and a y . We note that the ($) = ($*) at zero chemical potential. The 
numerical value of the pseudocritical temperature for various transitions in the QMS model and the PQMS model 
with the polynomial and the logarithmic potentials for the Polyakov loop has been given in Table [Vj It is evident 
from the table that the chiral transition gets shifted to the higher temperatures as a result of the inclusion of the 
Polyakov loop potential in the QMS model. 

We have chosen to compare the results of our calculation in the PQMS model with the corresponding results in the 
QMS model on a relative temperature scale T/T c x . Such a choice is justified on account of the Ginzburg-Landau 
effective theory, where absolute comparison of the characteristic temperatures between two models of the same 
universality class can not be made [3l|. The condensates start with fixed values a x = 92.4 MeV and a y = 94.5 
MeV at T = as shown in Fi^T] It is known from the lattice simulations that transition from hadronic matter to 
quark gluon plasma is a analytic and rapid crossover (0, The Polyakov loop potential inclusion in the QMS 

models makes the crossover in SUl{2) x SUr(2) sector quite sharp as the nonstrange condensate a x changes rapidly 
in the transition region. The Ua{X) anomaly does not cause any difference in the behavior of nonstrange condensate 
and a x remain unchanged in the presence as well as in the absence of V anomaly term. The variation of the 
strange condensate is lot more smooth on account of the large constituent mass of the strange quark m s = 433 MeV. 
The Polyakov loop potential inclusion has a strong effect on the strange condensate variation also and generates a 
significant melting of a y in our calculation. The interesting physical consequences of the earlier and significant melting 
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FIG. 2: The mass variations of the chiral partners as functions of reduced temperature (T/T* ) at z ero chemical potential 
(/i = 0), in the presence of axial Ua(1) breaking term, are plotted for (a, n) and (ao, r]') in Fig 2(a) and the corresponding 
mass variations, in the absence of the Ua(1) axial breaking term, are plotted in Fig 2(b) The dotted line plots are the mass 
variations in the pure QMS model, dashed- dotted line plots represent the PQMS:pol model results, and the solid line plots are 
the mass variations in the PQMSdog model. 



of the strange condensate will be an early emergence of mass degeneration trend in the masses of the chiral partners 
(K, k) and (77, /o) and an early setting up of a restoration trend on reduced temperature scale. In the presence 

of the Ua(X) anomaly, <j y temperature variation shows a little more decrease in the respective cases. 

Curves starting from the right end of the plot represent the variation of the Polyakov loop expectation value (3>) 
on the relative temperature scale at zero chemical potential. Though the thermodynamics of quark gluon plasma 
reproduced with the Polyakov loop polynomial potential is found to be in agreement with that of lattice simulations, 
upto twice of the critical temperature ([HI, O, HH) at higher temperatures (<£) increases above unity and this is 
unphysical. In the improved ansatz, logarithmic potential replaces the higher order terms of $ and <f>* in the 
polynomial potential by the logarithm of Jacobi determinant which results from integrating out six nondiagonal 
Lie algebra directions while keeping the two diagonal ones (l9l . l52| and thus the logarithmic divergence avoids an 
expectation value higher than one. This means that the logarithmic potential describes the dynamics of gluons more 
correctly and effectively. Keeping this in mind, we will be mainly focusing on the discussion of the results in our 
calculation with the inclusion of the logarithmic Polyakov loop potential, though the curves of the calculation with 
the polynomial Polyakov loop potential will also be shown. The real physical effect of the Polyakov loop potential 
inclusion in the QMS model on mcsonic excitations, will become apparent when the results of our calculation in the 
PQMS models are compared with the corresponding results in the QMS model. 



B. Meson Mass Variations 



We are calculating the masses of the scalar and pseudoscalar mesons at finite temperature in the presence of the 
Polyakov loop potential in the QMS model. We have collected the vacuum value of all the scalar and pseudoscalar 
meson masses in Table HU The mass modifications calculated at finite temperature (Eq |2"5| will be added to the 
vacuum masses of Table [TTJ The mass variations of the chiral partners as functions of reduced temperature, in the 
presen ce of axial Ua(^) breaking term, are plotted for (er, n) and (ao, r}') in Fig 2(a) and for (77, f n ) and {K, k) 
while the corresponding mass variations, in the absence of the !7a(1) axial breaking term, are plotted 
and Fig |3(b)[ Further, since the focus of our investigation is the influence of the Polyakov loop on the 



3(a) 



in Fig 
in Fig 

effective restoration of symmetries, we will be comparing the mesonic observables below and above 

In Fig 2(a){ the chiral partners (er, it) and (ao, rf) become mass degenerate in the close vicinity of reduced 
temperature T/T c x = 1. The masses of these particles are dominated by the contribution from the nonstrange 
quarks and a rapid crossover in the nonstrange sector (FigfTJ) appears as sharper and faster mass degeneration in 
our calculation in the PQMS model. Thus, the Polyakov loop inclusion in the QMS model makes a sharper mass 
degeneration as well as faster occurrence of chiral SUl(2) X SU r(2) symmetry restoration transition in the nonstrange 
sector. 

In Fig |3(a)[ the presence of the Polyakov loop potential in the QMS model generates, the similar trend of sharper 
and faster mass degeneration in the masses of the chiral partners (?7,/o) and (K, k). Though the mass degeneration 
of chiral partners (K, k) with 77 does not occur at T/T c x = 1, it sets up early in the PQMS models at T/T* = 1.3, 
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FIG. 3: The mass variations of the chiral partners as functions of reduced temperature (T/T x ) at zero chemical potential 
(/i = 0), in the presence of axial Ua{1) breaking term, are shown for (rj, fo) and (K, k) in Fig 3(a) and the corresponding 
mass variations, in the absence of the Ua (1) axial breaking term, are shown in Fig 3(b) The dotted line plots are the mass 
variations in the pure QMS model, dashed-dotted line plots represent the PQMS:pol model results, and the solid line plots are 
the mass variations in the PQMS:log model. 



while it occurs at T/T* = 1.5 in the QMS model. In the PQMS models, the intersection point of the fo and rj masses, 
occurs early when T/T x = 1.4, while in the QMS model it is found at T/T x = 1.7. This trend of mass degeneration 
reflects the effect of the Polyakov loop potential on chiral symmetry restoration in the strange sector and it results 
due to sharper and stronger melting of the strange condensate (FigfT]) in the influence of the Polyakov loop potential 
in the PQMS models. 

The Ua{]~) breaking generates the mass gap between the two sets of the chiral partners, (a, 7r) and (do, rj 1 ), i.e. 
Tot,- = m<7 < m aa = rrirji for T/T£ > 1. This mass gap results due to the opposite sign of the anomaly term (\Z2ca y ) 
in the scalar and pseudoscalar meson masses. Hence, it will be reduced due to the melting of the strange order 
parameter a y for the higher values of the reduced temperature T/T* > 1. Since the melting of the strange condensate 
is stronger and sharper (Fig|T]) in the PQMS model, the convergence in the masses of the two sets of chiral partners 
gets enhanced in these calculations. Thus the inclusion of the Polyakov loop potential in the QMS model also effects 
an early set up of J7a(1) restoration trend on the reduced temperature scale. 

Now we discuss the variations in the masses of chiral partners when the explicit Ua(^) symmetry breaking term 
has been taken as zero (c=0). We notice in Fig {2(b)1 and Fig {3(bjj again, the same sharper and faster trend of mass 
degeneration that we identify as the effect generated by the inclusion of the Polyakov loop potential. The rj meson 
degenerates with the pion in vacuum and stays the same for all temperatures in Fig 2(b) due to the absence of the 
anomaly term. Further the mass gap between the chiral partners (<r, 7r) and (ao, rj 1 ) becomes zero and all four of the 
mesons become degenerate at T/T x = 1.0. The T/T x numerical value, where the K, n, and rj masses degenerate in 
different models, is not influenced by the Ua{\) anomaly as expected since the nonstrangc condensate docs not have 
any anomaly dependence. Further, in Fig |3(b)| the intersection point of the fo and rj masses in the PQMS models, is 
obtained when T/T x is around 1.6 while in the QMS model, this intersection point is found around T/T* = 2.0. We 
are also obtaining the mild anomaly dependence of the intersection point of the fo and rj in all the models. Here, we 
note that the mass of fo in vacuum increases by about 60 MeV in the absence of anomaly. 

The temperature variations of meson masses, in general, result due to the interplay of the bosonic thermal 
contributions (decreasing the meson masses) and fermionic quark contributions (increasing the meson masses). Quark 
contributions which are negligible at small temperatures, dominate the mesonic contributions for high temperatures, 
and this generates a rising trend in meson masses, which ultimately leads to the mass degeneration of the chiral 



partners [38j. In the PQMS models, the one quark and two quark fermionic contributions are suppressed due to the 
presence of the Polyakov loop potential. Since the chiral phase transition is driven by the fermionic contributions, 
chiral restoration gets delayed due to the delay in the deconfinement transition and because of this we get higher 
value of the pseudocritical temperature T x in the Polyakov loop augmented quark meson linear sigma model. The 
higher value of the T x makes the ratio T/T x small. Hence, in comparison to the QMS model, the mass degeneration 
trend among the chiral partners, in general, sets up early in the PQMS models on the reduced temperature scale. 

The variation of meson masses with the polynomial Polyakov loop potential are similar though a little less sharp 
than the mass calculations with the logarithmic potential. The difference appears mainly because of difference in the 
Polyakov loop expectation value ($) with these two potentials. The calculations with the polynomial Polyakov loop 
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FIG. 4: The scalar 8s and pseudoscalar #p mixing angle variations with respect to the reduced temperature (T/T*) at zero 
chemical potential (fi — 0) are plotted. We have given the plots for the QMS, PQMS:log, and PQMS:pol models considering 
the cases in the presence as well as absence of the axial Ua(X) explicit symmetry breaking term. The dotted lines show the 
result with anomaly in the QMS model while the solid big dot lines show the result without anomaly. In the PQMS dog 
model, the thick solid lines represent the variations without anomaly while thin solid lines show the result with anomaly. The 
dashed-dotted lines are the variations with anomaly in the PQMS:pol model, while the dashed lines are the corresponding 
results in the absence of anomaly. 



potential make sense only for T < 2T*. 

The mass variation of scalar a and /o show kink around T/T c x = 1.8 in the QMS model while it is seen around 
T/T* = 1.4 in the PQMS models. The kink generation results because the meson masses seem to interchange their 
identities for higher values of the reduced temperature [38] . In order to have a proper perspective of the kink behavior 
in the curves, one has to study and analyze the scalar and pseudoscalar meson mixing angles. 



C. Meson Mixing Angle Variations 

The analysis of axial J7a(1) restoration pattern identification will become complete, only after studying the variation 
of scalar s and pseudoscalar Op mixing angles on the relative temperature scale in FigH] considering the cases in 
the presence as well as the absence of the axial J7a(1) explicit symmetry breaking term. The anomaly term has a 
strong effect in the pseudoscalar sector in the broken phase for T/T* < 1 while no effect of anomaly is found in the 
scalar sector. The nonstrange and strange quark mixing is strong, at T=0 one gets Op = —5°, which remains almost 
constant in the chiral broken phase. In the vicinity of T/T* = 1, the Op variations start approaching the ideal mixing 
angle Op — >• arctan-^ ~ 35° , the corresponding $p = 90°. Here, $p is the pseudoscalar mixing angle in the strange 

nonstrange basis (see Ref.[38[ for details). The smooth approach towards the ideal mixing in the QMS model, becomes 
sharper and faster in the PQMS model calculations due to the influence of the Polyakov loop potential. Further, the 
ideal mixing is achieved earlier on the reduced temperature scale in the PQMS models. In the absence of axial Ua{X) 
anomaly, the pseudoscalar mixing angle remains ideal Op = 35° everywhere on the reduced temperature scale. 

The 77 and 77' mesons become a purely strange rjs and nonstrange t]ns quark system as a consequence of the ideal 
pseudoscalar mixing, which gets fully achieved at higher values of the reduced temperature. In order to show this, 



we have plotted in Fig 5(a) the mass variations for the physical 77, 77' and the nonstrange-strange rjiys, Vs complex. 
Mass formulae m nNS and m vs are given in Table |TTJ Again, the smooth mass convergence trend, of the pure QMS 
model in — ^ m riNS and m ri — > approach, becomes sharper and faster around T/T* = 1 in the influence of the 
Polyakov loop potential in the QMS model. The exact m v i — > m rjNS and m v — > m vs mass convergence in the PQMS 
models, occurs closer to the value T/T c x = 1. 



13 



1100 



'Mill 



> 
o 

3 

CO 

S 700 



.-)(]() 



Ill 


// / / 
////,■■■'" 


S 

1J1S '< 




»7s '" 


V -^P 







0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 
T/T* 

■ QMS ■-■ PQMS:pol - PQMS:log 
(a)With Ua (1) anomaly term 





930 




920 


> 


910 






900 




890 








880 




S70 




860 



i i t i i i 




/ ' 

,/f PQMS:log 

T - - a 


j < 

J f 

/ * 
J 


1 — fo 

H &NS 

i I 

- i I 


QMS jfc 
iJl 
/• /* 

/* / 


r j 





1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 
T/T* 

(b)With U a(1) anomaly term. 



FIG. 5: Figure 5(a) shows the mass variations for the physical rj, rj' and the nonstrange-strange rjNS, rjs complex, on the 
reduced temperature scale (T/T c x ) at zero chemical potential (/i = 0). The masses of the physical a and /o anticross and the 
nonstrange-strange a ns — &s system masses cross in Figure 5(b) 



In Fig|4]for m a = 600 MeV at T = scalar mixing angle #s ~ 19.9° in the presence of anomaly, while 9s ~ 21.5° in 
the absence of anomaly. The 0$ around T/T* = 1 grows to its ideal value but for higher temperatures on the reduced 
temperature scale, in the chiral symmetric phase, the scalar mixing angle drops down to 9s ~ —51° and 9s ~ —54° 
in the respective cases considered with and without anomaly. In the presence of the ?7a(1) symmetry breaking term, 
this drop happens in the QMS model around T/T c x ~ 1.9 and due to the effect of the Polyakov loop potential the 
similar drop occurs earlier for T/T* ~ 1.5 in the PQMS model. In the close vicinity of these reduced temperatures, 
the m asses of the physical a and fo anticross and the nonstrange - strange (ctns ~ cs) system masses cross as shown 
in Fi ^5(b)| It means that after anticrossing the physical a becomes identical with the pure strange quark system 
as, while the physical fo becomes degenerate with the pure nonstrange quark system cttvs- A similar drop for the 
calculations without anomaly happens at a little higher value on the reduced temperature scale in respective models. 



VI. SUMMARY AND DISCUSSION 



We have calculated the meson masses and mixing angles for the scalar and pseudoscalar sector in the framework of 
the generalized 2+1 flavor PQMS model. We have used two different forms of the effective Polyakov loop potential for 
the calculation, namely, the polynomial potential and logarithmic potential. In order to investigate the influence of 
Polyakov loop potential on chiral symmetry restoration, these calculations have been compared with the corresponding 
results in the QMS model. 

The temperature dependence of nonstrange, strange condensates and the Polyakov loop field $ at zero chemical 
potential has been calculated from the gap equation in the QMS and PQMS models. Comparison of pseudocritical 
temperatures calculated from the inflection points of these order parameters indicates, that the chiral transition gets 
shifted to the higher temperatures as a result of the inclusion of the Polyakov loop in the QMS model. We further 
observe that the variation of the nonstrange condensate in the T/T* = 0.8 to 1.2 range becomes quite sharp due 
to the effect of the Polyakov loop potential in our calculation in PQMS models. We infer from the curves in the 
PQMS models that the inclusion of the Polyakov loop potential in the QMS model together with the presence of axial 
anomaly, triggers an early and significant melting of the strange condensate. The interesting physical consequences 
of the earlier melting of the strange condensate are an early emergence of mass degeneration trend in the masses of 
the chiral partners (K, k) and (rj, fo) and an early setting up of a Ua(X) restoration trend. 

The mass degeneration of chiral partners (a, tt) and (do, rj') in the close vicinity of T/T* = 1.0 becomes sharper 
and faster in our calculations in the PQMS model. This sharpening of the mass variations in the small neighborhood 
of T/T* = 1 results due to the stronger and sharper melting of the nonstrange condensate triggered by the presence 
of the Polyakov loop potential in the QMS model. Thus, we can corroborate also from the behavior of the chiral 
partners that the net effect of the Polyakov loop inclusion in the QMS model, is to make a sharper occurrence of the 
chiral SU(2)l x SU(2)ji symmetry restoration transition in the nonstrange sector. Further, the mass degeneration 
of chiral partners {K, k) with rj does not occur when the value of the reduced temperature is equal to one, it sets up 
early in the PQMS models at T/T* = 1.3, while it occurs at T/T* = 1.5 in the QMS model. In the PQMS models, 
the intersection point of the fo and rj masses, occurs early when the reduced temperature T/T c x = 1.3, while in the 
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pure QMS model this intersection point is found at T/T* = 1.7. This trend of mass degeneration emerges, again as 
a result of the sharper and stronger melting of the strange condensate in the influence of the Polyakov loop potential 
in the PQMS models. 

The J7a(1) breaking anomaly effect that leads to the mass gap between the two sets of the chiral partners, (a, it) 
and (ao, T)') i.e. m n = m a < m ao = m v > for T/Tq > 1, is proportional to the strange condensate a y . Since the 
melting of the strange condensate is stronger and sharper in the PQMS models, the convergence in the masses of the 
two sets of chiral partners will be enhanced in these calculations. Thus, the inclusion of the Polyakov loop potential 
in the PQMS models also effects an early set up of the Ua(1) restoration trend on the reduced temperature scale. 

The smooth approach of the pseudoscalar mixing angle Op towards the ideal mixing in the QMS model, becomes 
sharper and faster in the PQMS models due to the influence of the Polyakov loop potential. Further, in comparison 
to QMS model results, the ideal mixing on the reduced temperature scale is achieved earlier in the PQMS models. 
The 8s around T/T* = 1 grows to its ideal value but for higher temperatures on the reduced temperature scale, in 
the chirally symmetric phase, the scalar mixing angle drops down to 0$ ~ —51°. In the presence of Ua{X) symmetry 
breaking term, this drop happens in the QMS model for T/T* ~ 1.85 and in the PQMSdog model, the similar drop 
occurs for T/T c x ~ 1.5. In the close vicinity of these reduced temperatures, the masses of the physical a and /o 
anticross and the nonstrange-strange a^s — system masses cross. 
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